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Abstract 

A coarse-grained description of the restricted primitive model is considered in terms of the local 
charge- and number-density fields. Exact reduction to a one-field theory is derived, and exact 
expressions for the number-density correlation functions in terms of higher-order correlation func- 
tions for the charge-density are given. It is shown that in continuum space the singularity of the 
charge-density correlation function associated with short-wavelength charge-ordering disappears 
when charge-density fluctuations are included by following the Brazovskii approach. The related 
singularity of the individual Feynman diagrams contributing to the number-density correlation 
functions is cured when all the diagrams are segregated ito disjoint sets according to their topo- 
logical structure. By performing a resummation of all diagrams belonging to each set a regular 
expression represented by a secondary diagram is obtained. The secondary diagrams are again 
segregated into disjoint sets, and the series of all the secondary diagrams belonging to a given set 
is represented by a hyperdiagram. A one-to-one correspondence between the hyperdiagrams con- 
tributing to the number-density vertex functions, and diagrams contributing to the order-parameter 
vertex functions in a certain model system belonging to the Ising universality class is demonstrated. 
Corrections to scaling associated with irrelevant operators that are present in the model-system 
Hamiltonian, and other corrections specific to the RPM are also discussed. 
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I. INTRODUCTION 



The question of universality class of the critical point associated with a separation into 
ion-dilute and ion-dense phases in ionic systems has been a subject of an intensive debate 
for many years. The studies focused mainly on the restricted primitive model (RPM) and 
we shall focus on the RPM in this work. In the RPM equaly-sized hard spheres with posi- 
tive and negative charges of the same magnitute are dissolved in a structureless solvent (or 
in vacuum), and the system is globally charge-neutral. It was predicted theoretically |l| 
that the RPM can phase separate, and the separation was observed experimentally and in 
simulations. Stell put forward arguments that the RPM critical point belongs to the Ising 
universality class already many years ago 

sag 

, and the same conclusion followed later 
from the mesoscopic description described in Refs. 

urn 

Earlier simulation results for 

the RPM in continuum space were not sufficiently accurate to allow for definite conclusions 
concerning the universalityclass, however. Moreover, some early experimental results con- 
firmed the Ismg_criticality jsl, ll| , whereas other results were in conformity with the mean- field 



criticality 
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12j or indicated a crossover to the Ising-type critical exponents unusually 



close to the critical point 
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. In ternary aquaous solutions crossover to multicritical 



behavior was reported UM, and a tricritical, rather than a critical point was indeed observed 



in simulations 



eportea IJ 

Cm 



19, 20 



21 



22, 



23j in the RPM 



and predicted theoretically 
with the locations of the ions restricted to the sites of the simple-cubic lattice. The first 
firm prediction of tricriticaliy for that system appeared in Ref . [19j , where a strong argument 
based on symmetry considerations was given. 

The puzzling contradictions described above motivated careful theoretical, simulation 
and experimental studies of criticality in ionic sytems. In simple words, the Coulombic 
forces are long-range, and as such might lead to the mean-field (MF) criticality, as noted 



by Fisher 



24j . On the other hand, screeningeffects lead to short-range effective forces, and 



standard fluid criticality might be expected 



HQ 



One of the important issues is the range 



of the charge- density correlations in the critical region, i.e. beyond the infinite dilution 
regime J, 24 1. 

Subsequent precise experiments by Schroer, Wiegand and others [2 



31 



26 



27 



28, 



29, 



30 



showed that by careful identification and elimination of various additional effects 



(chemical instability, background scattering) Ising criticality is obtaind in all the systems 
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for which MF criticality or tricriticalty, or a delayed crossover were reported previously. The 
authors of Ref. repeated the experiments and found that the previous observations were 
associated with long- living met ast able states. After the equilibriation time (days) the Ising 



criticality was observed 



At present no RPM-like experimental ionic system deviates 



from the Ising criticality. Recent simulations also strongly support 

the Ising criticality. 

On the theoretical side the earlier arguments by Stell [3L M, and the field-theoretic de- 

n n 

termination of the universality class, briefly sketched in Refs. p, |6( are still considered 
as not fully satisfactory. Since exact analytical calculations are impossible, as is also the 
case for the three-dimensional Ising model, there is a demand for a theory analogous to 



the renormalization-group (RG) theory in simple fluids 4C| |4l( . The latter is restricted to 



short-range interactions, and is not directly applicable to the Coulombic forces 



In this work the theory introduced in Ref. 



5( is further developed, and exact derivations 



as well as technical details are explained on a formal basis. We start with the mesoscopic 



4^. In the 



description proposed in Reffl and described in more detail in Rets.Q 
coarse-grained description local distributions of the ionic species occur with a probability 
proportional to the Boltzmann factor exp(— /3f2 MF [0(x)p(x)]). Here f2 MF [0(x)p(x)] is the 
grand potential in the system with the charge- and number densities constrained to be 
0(x) and p(x) in an infinitesimal ('mesoscopic') volume cix centered at the space position 
x. It seems natural to reduce the two-field theory to a one-field theory depending only 
on the number- density of ions, which should be a proper order parameter for the gas- 
liquid type transition. However, in the RPM only charges interact, and it is not possible 
to perform such a reduction in a simple way. On the other hand, because there are no 
interactions between the masses in the RPM, the field p(x) can be integrated out easily. 
In sec. 3 we reduce the two-field theory to the one-field theory of the field 0(x), with the 
Boltzmann factor oc exp(— /37i e //[0]). The particular feature of the one- field theory is that 
the second functional derivative of T~C e ff[4>] assumes a minimum for the wavenumber k\, ^ 
in Fourier representation, and the higher-order terms oc J x A2n4 >2n {'x) have the property 
that for low valuse of the average number density A4 < and ^4 2 n > for n > 2. We find 
exact expressions for the number-density correlation functions in terms of the higher-order 
correlation functions for 0(x) in the one-field theory. 

In order to determine the scaling properties of the connected number-density correlation 



3 



functions we consider the corresponding Feynman diagrams in the (formal) perturbation 
expansion in *4 2 n- In sec. 4 we study the two-point charge- density correlation function G^. In 
the Gaussian approximation G^Ak) is singular for the wavenumber ^ 0. The singularity 
is removed beyond the Gaussian approximation by following the Brazovskii theory |43[ . We 
segregate all Feynman diagrams for any given connected number-density correlation function 
into disjoint sets according to a prescription described in sec. 5. By summing all diagrams in 
each set we obtain a finite expression, which can be represented by a skeleton diagram with 
lines representing the charge- density correlation function. We call such a series of diagrams 
'a secondary diagram'. In this way we cure the artificial singularity of individual diagrams 
following from the singularity of G^(kb). 

The secondary diagrams are again segregated into disjoint sets, according to their topo- 
logical properties, as explained in sec. 5. A series of the secondary diagrams belonging to 
a particular set can be represented by a hyperdiagram. The hyperdiagrams are associated 
with the same expressions for k — > as the corresponding diagrams in a model system with 
the Hamiltonian belonging to the Ising universality class, which contains also irrelevant 
operators. We could repeat the whole procedure of regularization of the hyperdiagrams, 
renormalization of the coupling constants and solve the flow equations describing the evo- 
lution of the renormalized coupling constants under changing the length scale, but this has 



already been done 



3, 41 1 for the model Hamiltonian yielding the same expressions for the 
vertex functions in the (formal) perturbation expansion. In sec. 5 we also describe corrections 
to scaling specific for the RPM. The last section contains a short summary. 

II. BACKGROUND 

A. Coarse-grained description for the RPM. 

In the field-theoretic, coarse-grained approach we consider local densities of the ionic 
species in mesoscopic regions dx around each point x, Pq,(x), where a = +, — . For fixed 
temperature and chemical potential the probability density that the local densities assume 



a particular form p+(x), p_(r) is given by 



p|>;(x),p* (x)] = H- 1 exp(-^[ P ;,p*_]), (1) 
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where 

S = / Dp* + J Dp*_exp(-f3Q MF [pl,p*_]), (2) 

(3 = 1/kT with k the Boltzmann constant and T temperature, and we introduced dimen- 
sionless densities 

Pa=PaO-+-, (3) 

where a + _ = (a + +a_)/2 is the sum of radii. In the above Vt MF [p* +1 p*_] is the grand potential 
in the system where the local concentrations of the two ionic species are constrained to be 
p* (r), pl(x), 

n MF [pi,p*_] = F h [pi, P *_] + u PM [pi,p*_]-p a [ P ;(x), (4) 

J X 

We use simplified notation J = J dx throghout the whole paper. In Eq.pjl is the 
Helmholtz free energy of the hard-core reference system, and we shall limit ourselves to the 
local-density approximation Fh[p* +1 p*_] = f d^fh(pK'x), pl(x)). In principle, more accurate 
approximations for -F/ l [p+, p*_] could be adopted. In the local density approximation 

' h(p) ' ( ) 

where the first term results from the ideal entropy of mixing, 

p* = p;+p*_, (6) 

and Ch{p*) is the volume integral over the hard-sphere contribution to the Ornstein-Zernike 
direct correlation function. 

In this work we focus on the RPM. In the RPM 5 2 Vt MF /8p* a 8p* p can be easily diagonal- 
ized due to the symmetry of the interaction potentials, and the eigenmodes have a natural 
physical interpretation. The first eigenmode is the number-density deviation of ionic species 
from the most probable value, 

ry(x)=p*(x)-p*, (7) 
and the second mode is the local charge- density in e-units (e = e+ = |e_|), 



(x)=p;(x)-p*(x). (8) 



The electrostatic energy is 



PU RPM {<P} = ^ I I 0( Xl )V(xi-x 2 )0(x 2 ), (9) 



B* = J_ = 

P T* Do' 



(10) 



V(x) = (11) 

and x = |x|. In the following, distances will be measured in cr + _ = o units. In Eq. (Q the 
contributions to the electrostatic energy coming from overlapping hard spheres are excluded 
because of the form of the potential V (see Eq. (jTT|) ). 

In the field theory introduced above the physical quantities are obtained by averaging 
over all fields p* +1 p*_ (or over <fi and r] in the RPM) with the Boltzmann factor and the 
grand potential Q is 

-/3fi = logS. (12) 

For small amplitudes of the fields and rj the functional (@J) can be expanded about its value 
Q^f F at the minimum, 

aq mf = n MF - n^ F = n 2 + n int . (13) 

Here fl 2 denotes the Gaussian part of the functional. In terms of the eigenmodes the 
Gaussian part of the functional (J 13)) assumes the form 



m2 = \ 



(14) 



where we simplify the notation for the integrals in the Fourier space f k = f j0]5- By using 
(0) we obtain 

Cl(k) = p* - 1 +(3*V(k), (15) 

and 



^(k) = 7o, 2 = /?|^U=o^=p 5 . (16) 
In the continuum-space RPM the Fourier transform of the potential (jlip has the form 



V(k) = (17) 



where k = |k| is in a -1 units. The wavenumber corresponding to the minimum of the 
potential V(k) is ^ 5]. The explicit form of 70.2 depends on the approximation for the 
hard-sphere reference-system; in any case, in the absence of short-range attractive forces it 
is a positive function of p^. Hence, rj is a noncritical field. On the other hand, ClAk) can 



vanish for k 0. The boundary of stability of AO in the (p*,T*) phase space is given 
by 

C^ih) = 0. (18) 

The critical fluctuations are thus 0(k&). 
The fli n t has the expansion 

/mw^i = £[ 7^^ 2m ^r- (19) 

where n denotes the summation with 2m + n > 3, and 72- m ,n denotes the appropriate 
derivative of f n at = and rj = 0. 

B. Ising Universality Class and Scaling 

The universality class of a given system is associated with a particular scaling form of 
the singular part of the thermodynamic potential in the critical region, or, equivalently, by 
the scaling properties of the large- distance part of the connected correlation functions. In 
particular, the Ising universality class is represented by the standard <p 4 theory with the 
Hamiltonian 

Hl[(f] = I l {t ° + fc2 M k M- k ) +J (- HoV(x) + | ^ 4 (x)). (20) 

The (renormalized) connected iV-point correlation function for the field if scales in the real- 
space representation according to [4^ 

G£( Xl , x,v; t, u, H) = t N f 3 G*(xit 1 ', ...,XNt u ; 1, u* , Ht~ A ), (21) 

where t, H and u are the renormalized coupling constants corresponding to the bare couplings 



the RG flow of u(£) upon rescaling 
4(] ] . In the case of a magnet t and H 



to, H and u respectively, and u* is the fixed point o 
the length unit, x — > x£ for £ — > 0, where u{£ = 1) = u 
correspond to reduced temperature and magnetic field respectively. The two independent 
critical exponents (3 and A are related to the other exponents via scaling relations [40] . In 
particular, /? = v{d - 2 + rj)/2 and A = u(d + 2 - 77) /2. 

Our purpose here is to determine the universality class of the RPM critical point. Thus, 
we have to study scaling properties of the correlation functions for the number-density field 
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r^(x) near the critical point of the phase separation, i.e. in the phase space region where 
((f)) =0. It is the correlation function for large separations, which may lead to divergent 
susceptibility and obeys scaling. In determining the universality class of a critical point one 
can separate the correlation function into a long- and a short-distance parts, and neglect 
the latter. In the next section we construct an effective one-field theory which allows for a 
determination of nonlocal parts of the number- density correlation functions. 



III. EFFECTIVE FIELD THEORY 



In this section we derive an exact effective Hamiltonian Ti e ff[4>], such that the grand 
potential is given by Q = —kTlog J D(j)exp(—/3T-C e ff[<j)}). Next we derive exact expressions 
for the nonlocal parts of the number-density correlation functions in terms of higher-order 
correlations for the field 0, with the Boltzmann factor exp(— /9"H e // [0])- We also introduce 
generating functionals for the relevant correlation and vertex functions. 



A. Derivation of the effective functional 

The functional AQ MF [<j),ri} (Eqs.JEiJ and ®) can be split into two parts, 

AQ MF [<f ) , V } = Q4<f ) } + Q r ,[<f ) ,ri} (22) 

where 

pnM = \(<p\cl^) + £ ^(i|0 2m ), (23) 

/ m=2 ^ m) - 

V] = 1(V\C°JV) + £ <0 2 "V>, (24) 

m,n ' 

and where Ylmn indicates a summation with n > 1 and 2m + n > 3. In the real space 
representation and r] are real functions. For convenience we use the Dirac brac-kets 
notation for the scalar product 

(f\g) = [ /*(x)y(x) = / /(k)fif(-k), (25) 

Jx Jk 

and \A\f) denotes the operator A acting on the state |/). In particular, when in real-space 
representation A has a functional form A(pc, x'), then 

|A|/) = £a(x,x')/(x'). (26) 



We introduce a functional of two external fields J(x) and v (x) by 

A~[J,v] = J DfaVWe-^+W' 



where 



v) = J Drje {vM e 



(v\ri) -/Xl^^r,} 



(27) 



(28) 



Because in the pure RPM (no other interactions included) C° (xi,x 2 ) = C° 5(xx,x 2 ) (see 
Eq.()16j)). Eq. (j24j) can be rewriten in the form 

(29) 



(30) 



Pttt,[<f>,V\= / /M,(0(x),7/(x)), 
Jx 

and we immediately obtain the exact result 



\ogE v [(f),v) = / w(0(x),«(x)) 



where 



u(<j>(x),v(x)) = log 



rf(r/(x))e- /3 ^ ( * (x) ' T?(x)) e 1,(x), ' (x) 



(31) 



The integrand u;(0(x),O) in Eq. (jSUjl can be expanded in a Taylor series with respect to 
0(x) and we obtain the formal expression 

/2m {Po) 



(0(x),o)=x; 



i2m/ 



(32) 



- (2m)! ^ * " 

where the explicit forms of /2m (Po) depend on the form of the Helmholtz free-energy density 
of the hard-sphere reference system. By substituting the above form of into Eq. (j27jl we 
obtain 

A~[J,0] = Z„H Q [J], So[J] = I D^e^ (33) 
where logZ v = J x /o(Po) (l°gAS[J, 0] is an extensive quantity), 



ICS 



m=2 



2m / ,2m, 



(2m)! 



and 



C° (x,x') + / 2 W(x-x' 



Am — 72m,0 + flm {Po}- 



(34) 



(35) 



(36) 
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FIG. 1: One-loop diagrams contributing to /2</> 2 (x). Thick lines and loops represent GjL and thin 
lines represent <p. Vertices with 2m thin lines and n thick lines represent 72m,n- 

Eqs. fj33|) - (}36|) define the effective functional. 

Eq.(|28j) can be written in a more general form, valid also beyond the RPM (i.e. for 
nonlocal CD, QQ 



Z® ■ exp 



l2m,n f j,2m/ \ ^ 



(2m)!n! J x 5w(x) n 
In the above the constant is Z® = J Dfe"^^ ^^, and 



e \{v\G° m \v)_ ^ 



Z v = Z° exp 



E 70,n 



70,n /" 5" 



n=3 



n\ J x (5w(x) 



eiH^I*)^. (38) 



Nonvanishing contributions to ^[0, 0] result from an even number of differentiations in 
Eq.flHH) with respect to v. In the case of the RPM (f\G° m \f) = G° m (f\f), and a double 
differentiation with respect to v gives a factor G° = 1/C° . The Wick theorem allows 
to represent 5^ [0,0] as a sum of all vacuum Feynman diagrams (no external points) with 
vertices 72m,n0 2m ( x ) such that n > 1 and 2m+n > 3, all at the same point x. From the vertex 
72m,n0 2m (x) there emanate n 77-lines, and all 77-lines have to be paired. Lines connecting 
different vertices as well as loops represent £r° . Finally, each diagram is integrated over x. 

log [</ >, 01 = J x cj(0(x), 0) is represented by all connected vacuum diagrams described 
above [4£j, |41[. co((j)(x), 0) is represented by the same diagrams, except that there is no 
integration over x. The term proportional to 2m (x) in the expansion of u;(0(x), 0) in Eq. (J32|) 
is represented by the sum of all connected diagrams, such that each diagram in the sum 
contains a certain number iV of vertices 7 2mj jnj .0 2m3 (x), and Ylj=i 771 j = m - I n particular, 
f2<f) 2 {x) is represented by a sum of all connected vacuum diagrams, where each diagram in the 
sum contains an arbitrary number of vertices 7o jn , and a single vertex 7 2) fc0 2 (x), where k is 
arbitrary. There are no zero-loop contributions to /20 2 ( x ) ; and the one-loop contributions 
are shown in Fig.l. The zero- loop diagrams contributing to A^^x) and As<f) e (x.) are 
shown in Fig. 2, and one loop diagrams contributing to ^.40 4 (x) are shown in Fig. 3. In 
these diagrams each thin line emanating from the vertex at x represents 0(x), thick lines 
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FIG. 2: Zero-loop diagrams contributing to ^4,40 4 (x) (a) and A§(ffi{x.) (b). Thick lines represent 
CrJL and thin lines represent <fi. Vertices with 2m thin lines and n thick lines represent j2m.n- 



FIG. 3: One- loop diagrams contributing to ^4</> 4 (x). Thick lines and loops represent GjL and thin 
lines represent <j>. Vertices with 2m thin lines and n thick lines represent 72m,n- 



represent G° , and the diagrams are multiplied by appropriate symmetry factors |40|, |41 |. 
A2m is represented by a sum of the corresponding diagrams with amputated 0(x)-lines. 

The above perturbation approach is convenient for a development of approximate theo- 
ries. Let us compare the above exact theory with the weighted-field (WF) approximation 



developed in Refs. 



44|. The WF approximation is of the mean- field type with respect to 
the field 77, because the effective functional is obtained by minimizing AQ MF with respect 
to 77 for a given field 0. TC e ff[(p] given in Eq. (j34j) has the same functional form as in the WF 



theory. Moreover, in the zero-loop approximation the term f% vanishes, and the coefficients 
Aim reduce to those obtained previously in Ref. Q|. Thus, in the zero-loop approximation 
the above exact effective theory reduces to the WF approximation, as expected. 
In the effective theory derived above we can rewrite S [J] in the form 



l [J] = Z^ l [Jl 



E 1 [J) 



, Z^m = 2 2m! Jx 



<5J(x)^ 



- e WK^\J) 



(39) 
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where = J Dfe sW^I/^ an d the charge-density correlation function in the Gaussian 
approximation, is inverse to C°^, 

G%(k) = l/C%{k). (40) 

By applying the Wick theorem we can represent Hi[0] by vacuum diagrams with 'hyperver- 
tices' from which there emanate 2m 0-lines. All lines are paired, and the line connecting 
vertices at xi and x 2 represents G^Jxi, x 2 ). The name 'hypervertex' used in Ref . comes 
from the form of *4 2m in the original, two-field theory (see Eq.fjHEJ) and Fig. 2), but in the 
effective one-field theory Aim are just usual vertices. For simplicity, in the rest of this paper 
we shall use the name 'vertex' for *4. 2m . 



B. Stability of the functional H e ff 

The effective functional of the form (jSljl was already studied in the WF theory j?| |44j , 
which is equivalent to the zero-loop approximation for / 2 and Am- Beyond the WF ap- 
proximation the Gaussian correlation function CrjL(k) diverges for the same wavevector 
k = k& 7^ as found in Ref. p , but along the spinodal line 

T* - -n* V ^ (411 

b ~ Po i + pMpIY (41) 

which is shifted compared to that found in the WF theory, where / 2 — (note that V^kjj,) < 
0). We assume that the denominator is positive for the relevant range of p$. For particular 
forms of fh this assumption remains to be verified. Divergent G°lAkh) for kf, ^ indicates 
that our theory belongs to the class of the Brazovskii field theories (431 ] . 

Another important property of our WF theory is the fact that for low densities, p$ < p^ c , 
the vertex A4 becomes negative, an d Am. > for n > 2 0, Q|. This property leads to a 



tricritical point on the sc lattice 



451 ] . and to a gas-liquid type instability in continuum 
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181 Ell- 



and some other lattice systems |45l. I46| . in agreement with simulation results 
Based on the WF and the simulation results we focus here on the general field theory 
laving the property that A^Pq) becomes negative for sufficiently low densities, with Aq > 
3, Q]- Within the exact theory developed above this assumption remains to be verified 
for particular forms of the reference system. Negative values of A4 for low densities bring 
our field theory in the corresponding part of the phase diagram beyond the class of the 
Brazovskii- type field theories, where A4 > 0. 
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C. Correlation functions in the effective field theory 



The correlation functions for the field can be obtained in a standard way from their 
generating functional AH [J, v] given in Eq. (J27j) . 

^ fiki+...+k n 



((^(xi)-/' l (x n ))) 



AH 5J( Xl ) fe i...5J(x T 



(42) 



and the connected correlation functions are given by 40j, |41 1 

£*!+••■+*» log AH 
= (JJ(xi) fc i...5J(x n )*"' J=0 ' ,,=0 " 



((0 fcl (xi)...0 fc "(x„))) c 



(43) 



((...)) denotes the average in the original two-field theory. From (j2*Tj) and (|3*3*|) we obtain 



((^( Xl )...0 fe "(x n ))) = (^(x 1 )...0*-(x n )>, 



(44) 



where (...) denotes averaging in the effective one- field theory. The connected correlation 
function (0 fcl (xi)...0 fc "(x ri )) con is represented by a sum of all connected diagrams containing 
external points Xi, ...x n from which there emanate hi, k n 0-lines respectively, and vertices 
.4.2m from which there emanate 2m 0-lines. All lines are paired. A line connecting a vertex 
or an external point at Xl with a vertex or an external point at x 2 represents GjL(x 1 ,x 2 ) 



40 



4l|. 



Let us consider the number- density correlation functions defined by 

1 <5 n AH 



AE5v(x. 1 )...8v(x n ) 



\j=0,v=0- 



From (}2Tj) we have 



where 



(( V {x 1 )...r l (x n ))) = i J D<Pe-^D n [^...^ n \ 



.D n [xi, -,x n | 



5v(xi)...5v(x. r 



\v=0- 



(45) 



(46) 



(47) 



By using Eq. (f3TJ) for local C° , in the case of xi ^ x 2 ^ ... ^ x n we obtain 



D n [xi,...,x n |0] = ^-GSexp 



E 



72m,n ^2m/ \ " 



(2m)!n! 



<fo(x) n 



u(Xi)...i;(x n ]ei«"'wN i ' l =0 . 



GS»H«>I 



(48) 

When Xj = Xj for at least one pair of external points, there are additional contributions 
to D n [xx, ...,x n |0] [40j, |4l|. We shall focus only on the case of x x ^ x 2 ^ ... ^ x n . From 
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the Wick theorem it follows that in the case of the RPM each diagram contributing to 
D n consists of n disjoint diagrams, each of them containing a single external point, and of 
vacuum diagrams. Thus, 

n 

A 1 [x 1 ,..,x n |0] = H 7? [0,O] nCW, (49) 

i=i 

where 

c(**) = £^ 2m w, (50) 

and a m are functions of p$ depending on the form of the free energy in the hard-sphere 
reference system. Each term a m 2m (xj) in Eq. ()5U|) is represented by the sum of all connected 
diagrams containing a single external point x», and an arbitrary number N of vertices 
72m J ,n J 2mj (xi) that satisfy ^2j=i m j = m - Eqs. (j4Hjl and (|5Uj) are of course immediately 
obtained by a direct differentiation of the exact formula (J3U)) . 

By substituting Eq.p9*j) into Eq.p5|l we obtain for xi 7^ x 2 7^ ... 7^ x n 

((7 ? (x 1 )...r/(x„))) = ±rf D0e-^//^C(xi)...C(x 7l ) = (C(xi)...CM>. (51) 
The connected correlation function for the number- density fields can be written in the form 



00 00 



((r ? (x 1 )...r / (x n ))) con =(C(xi)...C(x„)) co "= £ £ ^(0 2mi (xi)...0 2m "(x n )) co ". 

z — ' mi! z — ' m n ! 

mi=l m n =l 

(52) 

The above equation shows that for X] / x 2 / ... ^ x n the connected number-density 
correlation functions are given in terms of higher-order charge- density correlation functions 
with the Boltzmann factor oc exp f — /37i e ff[4>])- We introduce for convenience an additional 
field 

^(x) =a!0 2 (x), (53) 

where 

a! = <4(-2|i + ...), (54) 

and we rewrite Eq. fjB^j) as follows 



00 

(C(x 1 )...C(x n )) co " = £ ^ 7 ... £ ^ 7 ^ mi (xi)...^(x„)) co " = (55) 

* ' 1 777,1 * ' (7 , n T77.„ 



mi! ^ a x "m„. 

mi=l 1 m n =l A 

(^(x 1 )...^(x n )) con + ^(^ 2 (x 1 )...^(x n )) c - + 
af 
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Note that (£(xi)...£(x n )) con is given in terms of ( / 0(xi)...?/>(x m )) con with m > n, such that for 
m > n some of the external points are identical, and the number of distinct external points 
is n. (?/>(x 1 )...'?/>(x ri )) con is the dominant contribution to (((x.i)...({'X- n )) con in the phase-space 
region where the large-amplitude fluctuations are strongly damped by the Boltzmann factor. 

The phase separation, on which we focus in this work, is associated with a long-distance 
behavior of the number- density correlation function ((?7(xi)?7(x2))) cc "\ Contributions to 
((^(xi)?7(x2))) co?1 which are proportional to 5(xi — X2) are irrelevant when the behavior of 
((77(x 1 )^(x 2 ))) cc " 1 for |xx — x 2 | — > 00 is studied. Because of that we can approximate the 
exact correlation functions ((i](yii) . . .r](-x n ))) con by (C(xi)...£(x n )) con . By doing so we neglect 
only the local parts of the number-density correlations, because for xi 7^ x 2 7^ ... 7^ x n 
the above two functions are equal. Note, however that in the mesoscopic, or the coarse- 
grained theory we cannot describe the structure for distances comparable to the size of ions. 
Hence, the difference between ((77(x 1 )...r](x n ))) con and (((x 1 )...((x n )} con represents in fact 
the microscopic-scale contribution to the correlation functions, which is irrelevnt for the 
critical behavior. From now on we shall limit ourselves to the one field theory and to the 
correlation functions (V'(xi)...'^(x m )) con , which determine (C(xi)...£(x n )) con for n < to. 

D. Vertex functions and the free-energy functional in the effective theory 

In the one-field theory we introduce an effective generating functional for the connected 
correlation functions for the field ip that is defined in Eq.([53|h 



The functional ()56|) generates in fact also the connected correlation functions for the field (, 
according to Eq. (j55|) . It is convenient to introduce a functional analogous to the free-energy 
functional by a Legeandre transform of the above functional, 



where in Eq.()57|) V( x ) = ^^°E,'^eff[w}/Sw('x). The above functional can be expanded, 




(56) 



P T eff[lp] = logH e// H -(W\lp) 



(57) 



(3T eff [4j] = -(F 1 \^) + ^\C^) 
+ E n >J Xl --J x „^%^^(x 1 )..^(x„). 



(58) 
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is inverse to G^, where we introduced the notation 

G^(Ax) = (V(x)V(x + Ax)) co ™. (59) 

In diagrammatic expansion the vertex functions F n (x.i, x n ) are the one-particle irreducible 
(1PI) parts of the amputated n- point correlation functions, (V'(xi)...'^(x ri ))^ p , up to minus 
sign for n > 3. The 1PI diagrams cannot be split into two disjoint diagrams by cutting a 
single line. In the diagrams contributing to (^(xi)...-?/>(x n ))^ p the lines connected with the 
external points are amputated, 

^( Xl )..^(x n )) co " = / ... / ^(x')...^(x"))r p ^(x', Xl )...G^(x™,x n ). (60) 

The vertex functions F n (x!, x n ) determine the connected correlation functions 
(?/>(xi)...?/>(x n )) con , which in turn determine (C(xi)...('(x jn )) cc " 1 , i.e. the nonlocal parts of 
the connected correlation functions for the field i]. Thus, we have reduced the two-field 
theory to the effective one-field theory. 



E. Strategy in determining the universality class 

The universality class associated with the scaling properties of the correlation functions 
(if)(~Xi)...ijj(]L n )) con could be inferred from the form of the functional H\ip] such that 

log J Dije {wW e- mip] = logS e// H, (61) 

where logS e y/[w] is defined in Eq. (j56|) . Determination of the exact form of H\ip\ is not an 
easy task. Note, however that we are not interested in the exact form of the connected 
correlation functions, because scaling is obeyed only by their long-distance parts. The latter 
are generated by the singular part logS^j of the functional (J5T)Jl . 

logS e// H = logS'^H + logE^fw], (62) 

where log denotes the regular part of the functional. In a similar way the functional (JHEJ) 
can be separated into the singular and regular parts, T e ff[i/j] = rf^[^] + T^M], obtained 
by the Legeandre transform of logH^^ //M an d l°gSg^[u>] respectively. 

Our strategy in determining the universality class consists of the following steps. First 
the connected number-density correlation functions and their 1PI parts are found within 
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the perturbation expansion in the vertices Aim- Next we identify [■?/>], the singular part 
of T e ff[tp], and the vertex functions generated by rf^f^], i.e. those which are relevant in 
the critical region. In the last step we find a Hamiltonian H s [ip}, such that the Legeandre 
transform of the functional 

log J Dipe {wW e- msW (63) 

generates the same vertex functions. The scaling behavior of the vertex functions in the 
RPM is thus determined by the form of H. s [tJ)]. In sec. 5 we show that TC s [ip] has the form 
characterizing the Ising universality class. 

At the end of this section let us summarize the foundations of the effective theory. The 
nonlocal parts of the connected number-density correlation functions are given in Eqs. (j52j) - 
()55|) . The generating functional for (?/>(x 1 )...?/>(x n )) cc " 1 is defined in Eq. (j55j) . with the effective 
Hamiltonian Ti. e ff[<f\ given in Eq. (p£i|) . 



IV. ORIGIN OF THE CRITICAL SINGULARITY 



In the framework of the above coarse-grained approach the critical singularity of G^(0) 
was first found in Ref. in the WF approximation. The domain of validity of the equation 
for the spinodal line C^(0) = was not determined, however. Before describing the origin 
of the critical singularity, we shall briefly summarize the properties of the char ge-d ensity 
correlation function, G<^,(Ax) = (0(x)0(x + Ax)}™", determined already in Ref. 44j. The 
form of G^tf, is crucial for finding the domain of validity of the equation for the gas-liquid type 
spinodal line in the phase diagram. This is because the gas-liquid type phase separation can 
in principle be preempted by the charge ordering for G^(kfc) < 0, where C^ik) = 1/G<^(k), 
as is in fact the case on the simple-cubic lattice 3, 13 • 



A. in the self-consistent Hartree approximation 

The field theory with instability corresponding to kj, ^ was first developed by Brazovskii 
for the ip A theory. The Brazovskii theory leads to correct behavior of the correlation 
functions, i.e. the artificial singularity found on the Gaussian level is removed. For the RPM 
we keep terms up to O(0 6 ) in the functional (pMjl . because for low densities A± < 0, and 
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for stability reasons the 
with A G > 0. In Ref. 



unctional ()34j) can be truncated at the term oc 6 in the theory 
44j the Brazovskii theory for was considered for the RPM on 
the sc lattice, where a shift of the line of continuous transitions to the charge-ordered phase 
was found. Here we shall briefly describe the continuum case, where the above transition 
disappears. 

In the one-loop Hartree approximation the correlation function is given by an infinite 
series of effectively one-loop diagrams, shown in Fig.4 (top). In Fourier representation a 
single loop in the second diagram in Fig.4 (top) represents the integral 

Q = f G%(k), (64) 



k 



where the integrand G^(k) (see Eqs. (14 (J D and (J35j0 can be written in the form 

G° {k) = —~ . (65) 

In the above AV(k) = V(k) — V(kb), where for k = kf, V(k) assumes a minimum, and 



T = ^ ; / - + V(k b ). (66) 

Po 

Other diagrams in Fig.4 (top) are products of Qq. The self-consistent approximation is 
obtained, when in Eq. (JMj) the integrand is replaced by the correlation function which is the 
result of the whole resummation (Fig.4, bottom). The resulting equation is then solved self- 
consistently. In the above one-loop self-consistent Hartree approximation the k-dependence 
of the correlation function is the same as given in Eq. (j65|) . and only the critical parameter 
To is rescaled. We denote the rescaled critical parameter by r, and the correlation function 
in this approximation by G^. The corresponding k-integral of G^Jk) is 

e(T ),/ t ^ w = / t _|l_. (67) 



The self-consistent equation for G^{k) assumes the form (see Fig.4, bottom, and Ref. 

-l 



G H u {k) 



In the theory with the term oc 8 included in the effective Hamiltonian, there is another 
contribution, proportional to AsQ(t) 3 , in the above equation. From the definition (JBTIj) of r 
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FIG. 4: Top: a few diagrams contributing to the charge-charge correlation function in the one-loop 
Hartree approximation. Bottom: Diagrammatic representation of the self-consistent equation for 
the charge-charge correlation function. The thick line represents the correlation function G?±(k), 
thin line is the Gaussian correlation function (59, (&), the shaded circle and the bullet represent 
the vertices —A^ and —A§ respectively, and the k-integral is associated with each loop (see also 
Fig.2. 



we obtain the temperature corresponding to the singularity of G^Jkb) at r = (recall that 

v(h) < o), 

PoV(h) 



i + 0(0) f + #0(0) k (i + pS/ 2 (pS)) 



(68) 



The above temperature is larg er from zero for Q(0) assuming a finite value, for example in 



the RPM on the sc lattice 



4J. For 0(0) -> oo, however, T* — > 0. 



In continuum there are two sources of divergency of Q(t). The first, unphysical one is 
present for any r and comes from the behavior of V(k) for k — > oo. The charge- density waves 
with k — > oo would correspond to overlapping hard spheres, and should not be included. 



Following the standard procedure 



43, 



3, 



we expand AV(k) in a Taylor series about k = k b 
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and truncate the expansion in Ak = k — kj,. The resulting approximate form agrees with 
G^(k) for k fa kb, i.e. for dominant fluctuations, and its behavior for k — ► oo ensures finite 
values of Q{r). In the following we shall consider the integral (|67|) regularized as in Ref. 

3, 

s ^ = J - + ebs£ + 0(Ma) - ^/c^w) + (») 

where the subscript r stands for 'regularized'. Note that for reasonable regularizations (it 
is the neighborhood of kb that gives the relevant contribution to the integral) the value of 
Q r {j) depends on the regularization procedure only weakly (dots in the above equation). 

The other, physical divergency of Q r (r) occurs only for r = and is associated with 
the singular behavior of the integrand for k — > k^, i.e. with the dominant fluctuations 
0(k&) |44(. However, the temperature corresponding to r — > is T* — > (see Eqs.([5H)l 
and ([69)1 ). and for T* > we obtain in turn r > 0, hence the charge- density correlation 
function G^(k) remains regular for nonzero temperatures. Instead, a first-order transition 
is found for r < 

Q- In Ref.Q the fluctuation-induced first-order transition to the 
charge-ordered phase was identified with formation of an ionic crystal. The fact that the 
charge- density correlation function remains regular for nonzero temperatures is crucial for 
the occurrence of the gas-liquid type singularity in the continuum RPM. 



B. Critical singularity of 

In this section we describe the origin of the critical singularity of GW,(0), found already 
in Ref. 0. Moreover, based on the properties of G^ studied in Ref.|44| and summarized 
above, we show that the equation for the gas-liquid type spinodal line is well defined for any 
nonzero temperature. 

Let us first consider the correlation functions ( , 0(x 1 )...'0(x n )) con . The higher-order con- 
tributions to the number- density correlation functions in Eq. (J55J) will be considered later. 
Contributions to the n-point correlation function for the field ip are given by diagrams con- 
tributing to the 2n-point function for the field 0, with n pairs of external points identified, 
and multiplied by a", where a\ is given in Eq. flHH) . In particular, G^ (Eq. (|H9*j) ) is given by 
connected diagrams of the four-point function for the field 0, with the two pairs of external 
points identified with each other, and multiplied by a\. The contribution to G^(0) leading 
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to the critical singularity is given by an infinite series of diagrams which are of a form of 
chains of hyperloops connected by the vertices A4 5] (Fig. 5). The sum of the corresponding 
geometric series has the form 

G%(k) = A~g(k) [1 + A^k)}' 1 < (70) 

where g(k) is a Fourier transform of the function representing the hyperloop, i.e. a sum of 
all connected diagrams of the four-point function with two pairs of external points identified, 
which cannot be split into two distinct diagrams by splitting a single vertex A4 into two 
two-point vertices, and the sum is divided by 2 (symmetry factor). Another words, g(r), 
represented graphically by the hyperloop in Fig.5, has no contributions which are of a form 
of chains. We introduced a superscript b to distinguish the function given in Eq. ffTOj) from 
the exact form of the correlation function defined in Eq. (J59j) . 

Eq. (ffH|l stands a singular contribution to for k — ► either when g(0) — > 00 or when 

1 + .440(0) = 0. (71) 

The above can be satisfied for A4 < 0. Let us focus on the question for what temperatures 
g(0) is regular. The lowest-order approximation for g is represented by the first loop in Fig.5 
(bottom), and is given by 

~g (k) = I J drG^( r) V kr = \ J dk'G%(k')G%(\k - k'|). (72) 

G^Jk) is nonintegrable for r — > 0, because for r = G^Jk) diverges sufficiently fast 
when k — ► kf,. Thus, G Q A{k) is nonintegrable as well, and for r — » we have go(0) — > 
00. This singularity is associated with the charge- density waves with the wavenumber 
However, the divergency of G^ikb) is removed when the fluctuations are included within 
the Brazovskii approach described in the previous subsection. In a consistent approximation 
the hyperloop in Fig.5, bottom, can be written in the form 

g H (x!, x 2 ) = i(?^(xi, x 2 ) 2 + ...., (73) 

where G^(xi,x 2 ) is the charge- density correlation function, and dots indicate all remaining 
contributions. Because G^(k) is regular, g H (0) is regular as well, when the integrals are 
regularized for k — > 00 in the way described in the preceding subsection. Thus, the only 
singularity of G^(0) occurs at an infinite order in A4 and in this approximation is given 
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FIG. 5: Top: several diagrams contributing to the correlation function GjLr l (xi,X2), with the 
external points xi 7^ X2 shown as open circles. Black boxes represent a\. In the approximation 
corresponding to Ea. ((70|) . G^,(xi, X2) (represented by the pearl line) is given by an infinite series of 
diagrams which have a form of chains of hyperloops. Only diagrams with the hyperloops connected 
by the vertices A^ (shaded circle) are included in G^. Bottom: a few diagrams contributing to 
the hyperloop <?(xi,X2). The first contribution to <?(xi,X2) is go( x i> x 2) = G°, , (xi, X2) 2 /2. 

by Eq. ()71|) with g approximated by g H .The function inverse to that given in Eq. (j70|) can be 
expanded about k = 0, and we obtain 

C\k) = — L- = co + c 2 fc 2 + O(k'), (74) 

with Co and C2 expressed in terms of g H (0) and g H "(0) in a standard way. Note that the 
form of C b (k) for k — > is the same as the form of the inverse Gaussian correlation function 
in the Landau theory for simple fluids. The above result indicates that in the perturbation 
theory the separation into two uniform ion-dilute and ion-dense phases is present at an 
infinite order in A4 for Co = 0. The above discussion shows that in the continuum RPM this 
separation is not preempted by the charge ordering. 

We conclude this section by stressing that in the perturbation theory we should make 
a resummation of infinite series of diagrams twice. The first resummation, in the spirit of 



the Brazovskii theory |43j, cures the artificial divergency of the charge- density correlation 
function in the Gaussian approximation. Critical singularity is not exhibited by any indi- 
vidual diagram, but by an infinite series of diagrams having the form of chains of hyperloops 
(Fig. 5). The singularity is found only in the phase-space region where A\ < 0. Note that in 
the original functional AQ MF (Eq. ()13Jl ) the field 0(k&) is critical, and the field r] is noncrit- 
ical. Charge- density fluctuations lead to essentially different roles of both fields for A4 < 
- the field is turned to be noncritical, and the field r], or in fact ip oc <ft 2 , becomes critical. 
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FIG. 6: One-loop secondary diagram contributing to (0 2 (xi)0 2 (x2)) con (l.h.s.), and the correspond- 
ing series of Feynman diagrams (r.h.s). Thick and thin lines represent and G^, respectively, 
and open, gray and black circles represent external points, A4 and Aq respectively. 

V. UNIVERSALITY CLASS OF THE RPM CRITICAL POINT 

In this section we separate all diagrams contributing to the vertex functions for the field ip 
into disjoint sets according to the rules described below. The series of diagrams in each set is 
given by an expression that can be represented by a secondary diagram of a skeleton form. 
The secondary diagrams are again segregated into disjoint sets. The series of secondary 
diagrams in each set can be represented by a hyperdiagram. We show that the expressions 
representing the hyperdiagrams are the same as the corresponding expressions representing 
the diagrams in a certain model system belonging to the Ising universality class. The one- 
to-one correspondence between the hyperdiagrams representing the vertex functions in the 
RPM, and the diagrams representing the vertex functions in the model belonging to the 
Ising universality class indicates that the singular part of the free-energy functional has the 
Ising universality class form, up to additional terms associated with corrections to scaling. 



A. Secondary diagrams 

In the approximation studied in the previous section g H can be represented by a "sec- 
ondary diagram" of the same topological form as the first diagram on the r.h.s. in Fig.5, 
bottom, but with the thin line representing G^(xi, X2) replaced by a 'thick' line representing 
G r ^(x 1 ,x 2 ), as shown in Fig. 6. 

Similarly, all diagrams contributing to the connected correlation functions for the field 
V>(x) = ai</> 2 (x), can be divided into disjoint sets. Each set consists of one skeleton dia- 



gram 



4l| and of all diagrams obtained from this skeleton by insertions for any pair of points 



x',x", connected by a thin line representing GjL(x', x"), of subdiagrams which have a form 
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of diagrams contributing to the two-point function G^x.', x"). For each set a series of all 
its diagrams can be represented by a secondary diagram. The secondary diagram has the 
same topological form as the skeleton diagram, but all lines representing GjL are replaced 
by the lines representing G^. In a skeleton diagram there exists no pair of lines such that 
by cutting only these lines a subdiagram contributing to G^ could be extracted. Individual 
secondary diagrams are regular, because G^ik)™ is regular and integrable when G^ik) is 
regularized for k — > oo as described above. In the following we shall always consider sec- 
ondary diagrams. By doing so we automatically cure the divergency associated with the 
unphysical singularity of G^ikb). 

Because each vertex ^4. 2 m has an even number of 0-legs, and in diagrams contributing to 
(ip('Xi)...ijj('X2)) con an even number of lines representing G^ is connected with each external 
point, all diagrams contributing to (^(xi)...-?/>(x n )) con are 1PI. The secondary (skeleton) 
diagrams contributing to (?/>(x 1 )...'i/>(x2)) cori consist of loops containing n > 2 thick lines, 
and the loops are connected by the vertices A2m with m > 2. On the other hand, all 
diagrams contributing to (</>(xi)...0(x„)) con with xi 7^ ... 7^ x n are one particle reducible 
(1PR), i.e. two separate diagrams can be obtained by cutting a single line. 

In the following subsections we consider vertex functions and the free-energy functional 
in the perturbation expansion in the vertices A2 m - We shall segregate the corresponding 
secondary diagrams into different disjoint sets, associated with different contributions to the 
vertex functions. We focus on the general theory with A4 < and Ai n > for n > 3. We 
first analyze the lowest-order approximation. 

B. Perturbation expansion at the zeroth-order in the vertices Aim with m > 3 

Let us consider diagrams that contain loops connected only by the vertices A4. We shall 
first consider a general vertex part of the form of a loop with n > 3 lines. We show that such 
vertex parts are negligible in the critical region, and that the free energy-functional (J58j) has 
the same form as in a model system described by a certain Hamiltonian of the Gaussian 
form. 

Let us consider a general secondary diagram contributing to the n-th order vertex func- 
tion. At the one-loop order we have 

pone-loop^ _ ;Xn) K G ^ (xi _ X2 )... G ^( Xn _ Xl); (75) 
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n 

FIG. 7: A loop with n lines (n > 3) contributing to the n-point vertex function for 
the field if>, ^n( x ij ■■■■> x n)- Lines represent G<^. For the above contribution we have 
.L - Jx„ K°" e - /0 ° P (xi, ...,x n ) cx G„(0)« = 0. 

and the corresponding loop is shown in Fig. 7. For -F n (x 1; x n ) integrated over all arguments 
we find at the one-loop order 

/ - / Fr e -^(xi,-.,x n )oc 
/ ... / G^(x 1 -x 2 )...G^(x n -x 1 ) cxG^(0) n = 0. (76) 

The last equality follows from the fact that GjL(O) = (see Eqs.(JT5J) and (fPTj) ). and any 
diagram in the perturbation expansion of (7^(0) is proportional to some power of G^(0) = 
0. This is because all diagrams contributing to are 1PR. Note that the above property 
indicates short-range of the charge- density correlations in the uniform phase. For a secondary 
diagram with any number of loops the corresponding contribution to f ... j x F n (x.i, x n ) 
is also proportional to (5^(0) = for n > 3. 

Let us consider T e ff, the generating functional for the vertex functions F n (see Eq . (1581) ) . 
Because fluctuations relevant for the phase separation vary very slowly in space, let us focus 
on if>(r) = if) = const, in which case we have 

+ E„>J Xl --J x „^^r. (77) 

Let us consider the last term in the above equation. There are no zero-loop contributions 
to F n with n > 3 when the hypervertices A.% m with m > 3 are absent. As shown above, 
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at a higher number of loops the last term in Eq.(|77jl vanishes, and only the first two terms 
remain. From the resulting form of Eq. ()77|) it follows that neglecting the hypervertices A2n 
with n > 3 is analogous to the Gaussian approximation, provided that the fluctuations with 
wave numbers other than k = are disregarded. The role of the n-point vertex functions 
F n (ki, ...,k„), which for k { ^ behave as oc k\...k\ will be discussed in section IVEl where 
we briefly comment on all the irrelavant contributions to the vertex functions. 

Let us determine the contribution to coming from the secondary diagrams containing 
loops with more than two lines. As an example consider the secondary diagram shown in 
Fig. 8 (top). In Fig.8 (bottom) the hyperdiagram, denoted by /(xi,x 2 ), and obtained by a 
resummation of all diagrams of the above type is shown. We find 

/ / e t ^~^f(x 1 ,x 2 ) = f(k)V = (78) 

G^(k) 2 G H (k)G^(0) 2 f G^(x' 2 - x 3 )G^(x 2 - x' 3 ) = 0, 

^ x 2~ x 3 

where we used the fact that G^(0) = 0. From Eq. lfTKj) we see that the above contribution 
to G^(k) can be neglected. Similarly, all other diagrams of G^(k), containing loops with 
more than two lines, are proportional to some power of G<^(0) = 0. Hence, when the 
hypervertices A 2n with n > 2 are not included, the singular part of G^ is given by a sum of 
secondary diagrams of a form of a chain of loops with only two lines (Fig.9), and the inverse 
correlation function can be approximated by Eq. (J74|) . 

At the end let us consider the one-point function, i.e. the average density. The secondary 
diagrams contributing to (^(x)) = ai(0 2 (x)) at the zeroth order in the hypervertices A% 
and As are shown in Fig. 10, and the corresponding density shift is given by 

(V(x)> = ai g r (r) [l - A 4 ~g(0) + (-A^(O)) 2 + ...(-A 4 g(0))"] = Gj* (<>)*?, (79) 

where Eq.(|6T|) with G^L replaced by G </,</,, and Eq.(|70]) have been used, and we have intro- 
duced 

f ° m - f ° - Mi- m 

Note that F® (x) plays a role analogous to the external field. The above results lead to 



PreffM = -(FM + l^\C^\i;), (81) 
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FIG. 8: Top: A secondary diagram contributing to GU^ and containing the loop of the form (|75|). 
Lines represent G^, gray circles represent A4 and black boxes connecting the external points 
with the loops represent a\ ( Eg 154(1 . Bottom: A hyperdiagram obtained after a resummation of 
all diagrams of the form shown above. The corresponding expression in Fourier representation is 
given in Ea, (|78|) 



oocoococo = o^^o +c^[2) c CZ^ +C ^Z^ 3 CI) c CI^* :) + - 



FIG. 9: Secondary diagrams contributing to Gvw, (right-hand side) at the zeroth-order in the 
vertices A2 n with n > 2. GK ^ is obtained after a resummation of all diagrams of the above form, 
and is represented by the pearl line. Symbols are explained in the caption of the previous figure. 



FIG. 10: Diagrams contributing to the average density at the point shown as an open circle, and 
the hyperdiagram representing the sum of these diagrams. The black box represents a± fEq !54j) . 
The pearl line represents G^ (Eq.(J70J)) and Ff inside the box is given in Eq.(|HU|). 
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where F® and CvL(/c) are given in (|80p and (|74|) respectively. The same form of the free- 
energy functional is obtained in simple fluids with the coarse-grained Hamiltonian 

mW\ = lj c h w(k)${k)$(-k) - ^F{V(x) +(3H mt , (82) 

with the term 7i int neglected. 

C. Even and odd secondary diagrams 

The above observations indicate that the secondary diagrams contributing to the vertex 
functions for the field if) can be divided into two disjoint sets. In the first set, which we call 
the set of even diagrams, any contribution to the vertex function -F n (xi, x„) is represented 
by secondary diagrams such that any pair of points is either connected by an even number 
of lines, or is not directly connected. Such diagrams consist of chains of two-line loops, as 
in Fig. 9 and in Figs. 11-14 below. The diagrams in the other set contain loops such that 
there exist pairs of points connected by an odd number of lines. In particular, loops of the 
type shown in Fig. 7, where points x' and x" are connected by a single line representing 
G^(x' — x"), belong to this set of diagrams. We have shown that for if) independent of the 
space position such diagrams give a vanishing contribution to the free energy. In sec.5E we 
shall show that the odd diagrams in which some pairs of points are connected by three or 
more 0-lines lead to a modification of the coupling constants in the corresponding coarse- 
grained Hamiltonian (|K2*|) and to corrections to scaling, but do not change the universality 
class, which is determined by the even diagrams. 

Let us return to the function g representing the hyperloop (see Fig.5). Note that except 
from the secondary diagram of the same topological form as the first diagram shown in 
Fig.5, bottom, all remaining secondary diagrams contributing to g contain pairs of points 
connected by an odd number of lines and belong to the set of the odd diagrams. Hence, the 
approximation for g given in Eq. (f7H|) is justified as long as we limit ourseles to the set of 
even diagrams. 

The functional ()81|) was obtained by a resummation of all secondary diagrams with the 
two-line loops connected only by the vertices A±. The question which arises at this point is 
the role of the vertices Aq and As- For clarity of presentation we shall focus in the following 
subsection only on the even secondary diagrams. We shall determine the scaling behavior 
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of the corresponding contribution to the vertex functions. 



D. Even diagrams containing vertices Aq and As 

Let us first consider the three- and four-point vertex functions at the first order in the 
vertices Aq and As, and at an arbitrary order in A4. If we include only even diagrams, then 
the associated connected correlation functions are shown in Fig. 11 (left). By summing the 
infinite series of all diagrams of the above form we obtain the hyperdiagrams shown in Fig. 11 
(right). Any other secondary (i.e. of a skeleton form) diagram at the first order in Aq or 
.As, contributing to the respective correlation functions, belongs to the set of odd diagrams. 
The three- and four- point vertex functions at the first order in Aq and As respectively are 
just given by 

F°(k 1 ,k 2 ,^) = F° = ^ (83) 



and 



F°(k 1 ,k 2 ,k 3 M) = F° = ^. (84) 



Due to the translational invariance in the real space, the Fourier transforms are defined as 



n r r n 

(27r) d 5(^k i )F n (k 1 ,...,k n ) = / ... / F n (x 1 ,..,x n )exp(i^x J k i ). 



(85) 



The hyperdiagrams on the right in Fig. 11 are of the same form as the corresponding diagrams 
contributing to the connected correlation functions generated by the functional with 
the Hamiltonian (|82j). where 

mint = I (|W) + ^ 4 (r)j, (86) 

when calculated at the first-order in the couplings F® and F®. 

At the first order in As there exists a contribution to the connected three-point correlation 
function of the same form as shown in Fig.ll, top, left, but with the vertex Aq replaced 
by AsQrjl. However, the corresponding diagram belongs to the set of diagrams shown in 
Fig. 12, left, where the number of two-line loops representing g is arbitrary, from zero to 
infinity. Similarly, the contribution to the two-point connected correlation function with 
one vertex A4 replaced by AeQ r /2 belongs to the set of diagrams shown in Fig. 13, left, and 
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FIG. 11: Top: a secondary diagram contributing to the connected three-point correlation function 
at the first order in Aq (left) and the hyperdiagram resulting from the resummation of all diagrams 
of the above form. Bottom: a diagram contributing to the connected four-point correlation function 
at the first order in As (left) and the hyperdiagram resulting from the resummation of all diagrams 
of the above form. The black square represents As, the triangle represents the three-point vertex 
F3 and the open square represents the four-point vertex F°. The other symbols are the same as 
in the previous figures. 
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FIG. 12: Left: a secondary diagram contributing to the connected three-point correlation function 
at the first order in As and the hyperdiagram resulting from the resummation of all diagrams of 
the above form. The symbols are the same as in the previous figures. 

the corresponding contribution with one vertex Aa replaced by A^Qljl? belongs to the set 
of diagrams shown in Fig. 14, left, bottom. 

Let us consider an arbitrary vertex function for the field ip, and focus on a particular 
contribution of a form of an even 1PI secondary diagram containing a given number of the 
vertices Aq and As- Such a diagram consists of these vertices connected by chains of two-line 
loops, and the loops in the chain are connected by Aa. Some of the chains of loops may 
end with the loop representing Q, as in Fig. 12, left. By changing the number of loops in 
the chains we obtain other diagrams. All diagrams obtained in this way form a set. By a 
resummation of all diagrams belonging to this set we obtain a contribution to the considered 
vertex function, which can be represented by a hyperdiagram. The topological form of the 
hyperdiagram is given by any member of the set, if the chains of loops connecting the vertices 
Aq and/or A% are replaced by the pearl lines, as shown in Figs. 11-15. 

It is instructive to consider a few examples. The even secondary diagrams contributing 
to and containing a single vertex Aq and As, are shown in Figs. 13 and 14 (left), 
respectively. By summing the infinite series of all diagrams of the above form we obtain the 
hyperdiagrams shown in Figs. 13 and 14 (right), with the pearl lines representing the 'bare' 
correlation function 1)70)1 . The diagrams contributing to the connected four-point function 
at the second order in As are shown in Fig. 15. 

All even secondary diagrams can be separated into disjoint sets obtained in the way 
described above. By a resummation of all diagrams belonging to any given set we obtain a 
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FIG. 13: A diagram contributing to the number-density correlation function at the first order in 
^.6 represented by the bullet (left). The contribution to G^, given by the infinite series of such 
diagrams with all numbers of hyperloops, is shown schematically (right) as a hyperdiagram, where 
the pearl line represents the bare function (|7U[) and the open triangle represents — F® (see Fi dIVB| . 

hyperdiagram. For any hyperdiagram contributing to the vertex function in the RPM there 
exists a diagram of the same topological structure and representing the same expression, 
contributing to the vertex function in the model given by the Hamiltonian (|H2*|) and (|86|h 
and vice versa. This is just a topological property related to the possibility of connecting 
three- and four point vertices by lines, when the functional forms of the pearl line in the 
hyperdiagrams and the line in the diagrams corresponding to the above Hamiltonian are the 
same. We have thus reduced the question of the critical behavior of the vertex functions 
determined by the even diagrams to the question of the critical behavior of the system 
described by the Hamiltonian (JH2J) with (jSo^l. where C^(k), F°, F° and F® are given in 
Eqs. (|71j). (JHUj) . (jHSj) and (|83jl respectively. We could include also higher-order vertices 
•4.2m with m > 4, and we would obtain in the same way additional terms oc ip m m Eq.( 
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However, such terms are irrelevant in the RG sense 
By shifting the field, 

we can remove the term oc ip 3 in Eq. 1)86)1 . and we obtain the standard form of the coarse- 
grained Hamiltonian representing the Ising universality-class, Ti[ip\/c2 = 7~i- I \p\ + const., 
with 'H I [^p\ given in Eq. (J20j) . where the two relevant scaling fields are 

/ ^0 ^03 \ 
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FIG. 14: Left: secondary diagrams contributing to the number-density correlation function at the 
first order in As, represented by the black square. The contributions to G,^, given by the infinite 
series of such diagrams with all numbers of loops representing g H , are shown schematically (right) 
as hyperdiagrams, where the pearl line represents the bare function (J7UJ) and the open square 
represents — F4). 

and 

to= ( c o~^> k 1 , (89) 

and Mo = F4C2 1 . Note that the scaling fields H and to are related to the temperature and 
the chemical potential (through p$) in a somewhat complex way. Moreover, the values of 
the above parameters depend on Q r and g(k), which in turn depend on the regularization 
procedure of the k-integrals of and its powers. We shall not study the nonuniversal 
properties in this work. 

The contributions of all even diagrams to the vertex functions for the field ip are of the 
same form as the vertex functions determined by the Hamiltonian given in Eq. (j2()j) and 
belonging to the Ising universality class. Hence, the scaling properties are also the same. 
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FIG. 15: Top: A secondary diagram contributing to the connected four-point number-density 
correlation function at the second order in ,4s (black square). Bottom: hyperdiagram representing 
an infinite series of all such diagrams. The pearl line represents the bare function l|70|l. 

We need to find out whether the odd-diagrams contributions to the vertex functions can 
alter the scaling properties. 

E. Odd diagrams containing vertices Ae and As 

In this section we consider the secondary diagrams contributing to the vertex functions 
for the field ip oc 2 , in which there exist pairs of points connected by an odd number of 
lines representing G^. We focus on the theory with the vertices A\ < and ^4. 2n > for 
n > 3. Our purpose is to show that the odd diagrams give the contributions to the vertex 
functions which either scale in the same way as the even diagrams, or are associated with 
corrections to scaling. Because an even number of the 0-lines emanates from each vertex 
and each external point, and an even number of lines is amputated from the external points 
in diagrams contributing to the vertex functions for if), the vertices connected by an odd 
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number of lines form closed loops. 



1. Irrelavant odd diagrams 

In sec.5.B we have already considered loops formed by single lines, as shown in Fig.7, 
where F„(ki, ...k n ) oc Yli G^ih) oc kf...kf...k^ for ki — > 0. Let us focus on a contribution 
to F n (ki, ...k n ) that is of the form F^(ki, ...k n ) oc k 2jl ...k^ ji ...k%> n , where ji = 0,1 and 
Y2i ji > 0, and in the corresponding secondary diagrams no pair of points is connected only 
by a chain of the two-line g-loops, and no loops representing Q are present. A series of such 
secondary diagrams for a given {ji} can be represented by a hyperdiagram containing no 
pearl lines. We have J ... J x F?^(xi, ...x n ) = 0, because ji > 0. 

By following the considerations described in the preceding subsection, we can segregate 
all secondary diagrams contributing to a particular vertex function into disjoint sets. In the 
secondary diagrams belonging to each set the vertices As, As and F?^ are connected by 
the chains of the two-line g-loops. All diagrams in any given set are of the same topological 
structure, when the chains of the two-line loops are replaced by the pearl lines. The series 
of all diagrams in the given set is represented by a hyperdiagram. There is a one-to-one 
correspondence between the hyperdiagrams described above, and diagrams obtained in the 
perturbation expansion in the theory with the coarse-grained Hamiltonian (|82|) and with 7ii n t 
(|Kfi|) supplemented with contributions of the form 0(ip) oc J k ... f k 5(^2™ kj) Yl™ tp(k.i)k^ 1 
with n > 3. Note, however that in the critical region the terms 0(ip) correspond to irrelevant 



operators in the RG sense 



4l|. 



2. Relevant odd secondary diagrams 

Let us consider odd secondary diagrams contributing to the vertex functions for the field 
t/>, other than those proportional to G^(0) = 0, which are irrelevant. In such diagrams 
there exist closed loops in which neighboring pairs of points, xi and X2, are connected by 
subdiagrams contributing to G^xi — x 2 ) = (0 3 (xi)0 3 (x 2 )). The diagrams contributing to 
G^ipci — x 2 ) are shown in Fig. 16. Because G^ik) ~ k 2 for k — > 0, and the integrals are 
regularized for k — >• oo as described in sec. 4, the individual diagrams contributing to G^{Q) 
are finite. The series of diagrams of the form shown in the second diagram on the r.h.s. in 
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FIG. 16: Relevant odd secondary diagrams contributing to the correlation function Gs(xi,X2) = 
(^> 3 (xi)(/> 3 (x2)}. The dashed line represents the series of all diagrams of the form shown on the 
r.h.s. 

coo of x> °o~ >o co 



FIG. 17: Top: Relevant odd secondary diagram contributing to the two-point correlation function 
for the field ip. The dashed lines represent the sum over all sub diagrams contributing to G3, and 
the solid lines represent G^. The gray circles, black squares, open circles and black boxes represent 
Aa, As, the external points, and a\ respectively. Bottom: secondary diagrams contributing to the 
hyperloop representing g, when the even as well as the relevant odd diagrams are included, and 
the vertices A2m with m > 4 are neglected. 

Fig. 15 is also regular for the same reason. In the theory with A% > the series of all diagrams 
shown on the r.h.s. in Fig. 15 is regular for k — * 0, and we find G 3 (k) = g® + g\k 2 + .... Note 
the significant difference between the pearl line (Fig.9) obtained after the resummation of 
the even diagrams, and the dashed line obtained after the resummation of the odd diagrams 
(Fig. 16). Because A4 < and Aq > 0, only the pearl line can be singular for k — > 0. 
In Fig. 17, top, we show a relevant odd secondary diagram contributing to the two-point 
correlation function for the field ip. This diagram represents in fact the series of all secondary 
diagrams that contain subdiagrams contributing to G 3 (x — x') for each pair of the vertices 
x and x' connected by the dashed line in Fig. 17. We shall keep the name 'a secondary 
diagram' for such a series of secondary diagrams. The corresponding secondary diagrams 
contributing to the hyperloop g (see sec.IVB, Fig.5) are shown in Fig. 17, bottom, where the 
irrelevant diagrams are not included. In this approximation we obtain for G b ^{k) the same 
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FIG. 18: Relevant odd secondary diagram contributing to the re-point vertex function for the field 
ip. The black squares represent As- The dashed line represents the three-point correlation function 
G3. Thin lines represent the (amputated) 0-lines. 

expression as given in Eq. ffTO]) . but with g(k) given by 



where g H is given in Eq. (j73|) . and G3 is shown in Fig. 16. For k — > we obtain the bare two- 
point correlation function G^Ak) of the same form as given in Eq. (j74J) . but with modified 
coefficients cq and C2. Note that the secondary diagrams contributing to the bare function 
GcLAk) contain vertices other than A4, unlike in the case of the even diagrams. Note also 
that the secondary diagrams which contain chains of loops such that the two ends of the 
chain emerge from the same vertex, or one end of the chain represents Q{r) (as is the case 
in Fig. 14), are not included in the bare function G b ^. As in the case of the even secondary 
diagrams, only linear chains, with no chain-like branches, contribute to G^. 

In Fig. 18 we show the contribution to the re-point vertex function for the field ip, 
F° dd (k!,...,k n ), of the form Ui^ih). Note that G 3 (k) = g° + g\k 2 + 0(k 4 ) is finite 
for k —>■ 0, therefore the above contribution to F° dd (0, 0), being of the form oc G^(0) n , 
is nonvanishing as well. Hence, the above contribution to F n is relevant in the critical re- 
gion, and should be taken into account. By F® (ki, k„) = F® odd + 0(k 2 ) we denote the 
contribution to F° dd (ki, ...,k„) which is given by all odd secondary diagrams such that no 
pair of vertices is connected only by a chain of the g-hyperloops and no loops representing 
Q are present. The vertex part F® odd (0, 0) = F® odd is relevant in the critical region, and 
pOodd gho^Qo! Iqq taken into account in the same way as the vertices Am are. Thus, the 
contributions to any vertex function, which are relevant in the critical region, consist of the 
even secondary diagrams, as well as of the relevant odd secondary diagrams (those which 




(90) 
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contain the dashed lines). Any relevant 1PI secondary diagram contains vertices Aq, As 
and the vertex parts F® , F® odd , from which there emanate 3 or 4 linear chains of hyper- 
loops respectively, and the chains are of the form shown in Fig. 17. An example of an odd 
secondary diagram contributing to the four-point connected correlation function is shown in 
Fig.19. 

As in the case of the even diagrams, we can segregate all the relevant diagrams (both 
even and odd) into disjoint sets. To a given set there belongs a diagram in which particular 
pairs of vertices or the vertex parts are connected by the linear chains of the hyperloops. 
All secondary diagrams obtained from this particular diagram by changing the number of 
the hyperloops in the chains belong to the same set. This is analogous to the segregation 
of the even diagrams described in sec.VD. There is one-to-one correspondence between the 
vertex functions in our theory and the vertex functions obtained in the theory with the 
coarse-grained Hamiltonian (|82jh when the terms 
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3 3 

i / ... / (2n) d 5(^k l )F^ dd (k u ...M)flm 



and 



1 

4!./ k 



I (2n) d 6(J2 k^r d (ki, k 4 ) J] ^(k,) 

k 4 



(92) 



are added to Hint (Eq. 
relevant contributions, 



Because of the form of F° odd (ki 



■j k n ), 



we can separate the 



pQodd 

~1T 



3 3 pOodd 

(2n) d 5(Y,K)Hm) ~ 



3! 



V>(x) ; 



(93) 



and 



pOodd 



ki 



k 4 



4 4 pOodd 



4! 



^(x) 



(94) 



In the remaining contributions the integrands are proportional to 0(kf) for ki — > 0. Such 



contributions to Timt are irrelevant in the RG sense |4l|. The vertex functions for the field 
if) scale in the same way as the vertex functions in the model given by the Hamiltonian 
(I82|) and (J86|) . but with all the parameters modified according to the above discussion. In 
particular, F° -> A n ja n x + F® odd . 

With higher-order vertices Aim, rn > 4 included, similar vertex parts with loops made 
of subdiagrams contributing to (</> 5 (x^^ 5 (x 2 )) con etc. would contribute to F 3 °, F° in Eq. 
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FIG. 19: Top: A relevant odd secondary diagram contributing to the four-point connected corre- 
lation function for the field ip. The gray circles, black squares, open circles and black rectangles 
represent Aq, As, external points, and a\ respectively. The solid and dashed lines represent the 
two- and the three-point correlation functions for the field (f), and G3 respectively. All diagrams 
obtained by changing the number of the solid- and the dashed-line loops in the four linear chains 
belong to the same set of diagrams. Bottom: A hyperdiagram obtained after the resummation 
over all numbers n > 1 of the solid-line loops and over < m < n — 1 dashed-line loops in each 
linear chain, and over all subdiagrams contributing to the bare four-point vertex functions. The 
pearl-line represents the series of all linear chains and the bare four-point vertex functions are 
represented here as the shaded squares. See text for more datails. 



and to the irrelevant operators. Modification of the coupling constants affects only 
the nonuniversal properties of the theory. Corrections to scaling associated with irrelevant 



operators have been studied before |4lj, and are not specific for the RPM. 

Let us summarize the above results. From all the secondary diagrams contributing to the 
vertex function for the field ip, F n (ki, k n ), we can distinguish a set of secondary diagrams 
contributing to the 'bare' n-point vertex function. To the above set there belong the vertex 
-4-2rt/(Xi and the odd secondary diagrams such that no pair of vertices Aim is connected by the 
linear chain of the hyperloops representing g (including the one-loop chains representing g H ). 
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Also, no loops representing Q are present in the diagrams contributing to the bare vertex 
function. The bare vertex function has the form i^ are (ki, k n ) = Ami cii+F® odd +0(k 2 ) for 
k -> 0, where F® odd is a real number. The remaining secondary diagrams contributing to the 
vertex function F m (ki, ...,k m ) consist of the subdiagrams of the form described above, and 
these subdiagrams are connected by the chains of the hyperloops representing g (including 
the one-loop chains representing g ). All the secondary diagrams contributing to the n- 
point vertex function for the field ip can be segregated into disjoint sets. Each set contains 
all the secondary diagrams obtained from one particular secondary diagram by changing 
the number of the hyperloops in the chains connecting the subdiagrams that belong to the 
set representing the bare vertex functions. Also the secondary diagrams obtained from the 
chosen diagram by changing the subdiagrams that contribute to the bare vertex functions 
belong to the considered set. All the secondary diagrams in any given set are of the same 
topological structure when the linear chains are replaced by lines, and the subdiagrams 
contributing to the bare n-point vertex function are replaced by vertices from which n lines 
emanate. The series of all the secondary diagrams in a given set can be represented by 
a hyperdiagram. The hyperdiagram has the same topological structure as the secondary 
diagrams in the set. In the hyperdiagrams the functional form of the line representing 
the sum over all numbers of the hyperloops in the linear chain is l/(co + 02k 2 ) in Fourier 
representation, and the n-point hypervertices are of the form A2nl^\ + F® odd + 0(k 2 ) for 
k — > 0. The generating functional for such vertex functions is obtained from the Hamiltonian 
f!82jl with (jHHJ), that contains also additional, irrelevant contributions to Hint that yield 
corrections to scaling. 

F. Corrections to scaling specific for the RPM 

We found that the connected correlation functions for the field ip scale in the same way 
as in the Ising universality class. However, the correlation functions for the number-density, 
given in Eq.(|55|). contain also other contributions. In particular, the one- and two-point 
number-density correlation functions in the RPM assume the forms 




(95) 
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FIG. 20: Left: Secondary diagrams contributing to the next-todeading order term in the expansion 
of (C( r )} fEa. (|95[) ) at the zeroth order in the hypervertices A^n with n > 2. The external point is 
shown as an open circle, gray boxes represent 02/2, lines represent and gray circles represent 
A4. Right: Hyperdiagrams resulting from the resummation of all diagrams of the type shown on 
the left. 

and 

(,(r 1 Mr 4 rH#i)^)r+ E E »STt ^ wn^r" (96) 

, „ Oi n\m\ 

n,m>l n+m>2 1 

respectively, where the averages are obtained with the Hamiltonian (|H4j) . In Figs. 20 and 
21 we show the secondary diagrams contributing to the next-to- leading order terms in the 
above expansions. On the right the hyperdiagrams obtained after the resumation of all such 
diagrams are shown. In general, the leading-order correction to the iV-point connected 
correlation function for the field ip is given by the iV + l-point connected correlation function, 
where two of the external points are identified. Scaling forms of the connected correlation 
functions, Eq. (J2T|) . give 

G#( ri , r N ; t, u, H) = \t\ Np G Na (T X /Z, r N /& H\t\-^) + (97) 
|t| (Ar+1) ^ +1)s ( ri /e, n/e, r*/fc H\t\- A ) + 0{\t\W), 

where Gn s (t 1 /^, r N /^; H\t\~ A ) is the scaling function. Consider now the correction term 
to the grand potential. Because the relative correction to any connected correlation function 
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FIG. 21: Left: secondary diagrams contributing to the next-to- leading order term in the expansion 
of G^ (Eg .([96))) at the zeroth order in the hypervertices A n with n > 4. The external points are 
shown as open circles, black and gray boxes represent a\ and 02/2 repsectively, lines represent 
Gsch and gray circles represent A4. Right: Hyper diagrams resulting from the resummation of all 
diagrams of the type shown on the left. 



is oc {t^, the relative correction to their generating functional has the same scaling form, 
and we obtain 



Q(t,H) = \t\' 2 - a iu s (H\t\- A ) + \t\^uo cs {H\t\- A ) + ... 



1/3, 



All derivatives of the grand potential can be written in the form 



X oc \t\ a 



i + |t|^ M Off|t|- A ) + o(|t| 



2/3 A 



(98) 



(99) 



where a is the corresponding critical parameter in the Ising universality class. The leading- 
order correction to scaling, specific for the RPM, is given by the exponent /3. 



VI. SUMMARY 



The theory outlined above is relatively complex compared to the theory of critical phe- 
nomena in simple systems. This complexity follows from the indirect nature of criticality 
in ionic systems. The long-range critical number- density fluctuations in charged systems 
are not induced by mechanical interactions. Unlike in systems with short-range interac- 
tions, the phase separation is induced by strong charge- density correlations. Coulombic 
forces support formation of charge-ordered structures with a microscopic distance between 
oppositely charged neighbors. Individual microscopic states are charge-ordered with a high 
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probability as predicted theoretically j^J and seen in snapshots 5^, 51, 0], but fluctu- 
ations restore the uniform structure at larger length scales (and large observation times). 
Mathematically the restored disorder is described within the Brazovskii approach. In our 
case we performed a resummation of singular Feynman diagrams (singularity resulting from 
charge-ordering) to obtain regular secondary diagrams. The charge-ordered 'living' clusters 
that are formed in different microstates interact with each other with short-range forces, 
and this observation suggests standard criticality. In microscopic description interactions 
between clusters of various sizes, shapes and orientations should be considered. The above 
complexity in our theory is reflected in the absence of critical singularity in individual sec- 
ondary diagrams. The critical instability of the uniform phase occurs only when an infinite 
series of secondary diagrams contributing to the correlation function (0 2 (x 1 )0 2 (x 2 )) con is 
calculated in the perturbation expansion. In the diagrams of a form of chains of n loops the 
correlations with n — 1 intermediate points are included. The series of chains of loops plays 
an analogous role as the Gaussian correlation function in simple fluids. The hyperdiagrams 
with the lines representing such series play in turn an analogous role as diagrams in the 
standard theory of critical phenomena. 

From the point of view of the field theory, the above work concerns the theory with the 
action of the form ()34|) with Aa < and A% n > for n > 3, and in Fourier space the action 
becomes unstable for the wavenumber k b ^ 0. The key properties A$ < and Am > have 
been observed within the WF approximation at low densities for the RPM with the ideal 
entropy of mixing and with the Percus-Yevick [3] reference system. It is plausible that 
for the exact form of A± the condition Aa < is satisfied in the RPM for low densities, but 
it remains to be proven. This work shows that the model for which the above is satisfied 
belongs to the Ising universality class. 
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